﻿/*******************************************************************************
 * Copyright (c) 2007
 * Sergey Berezin, Dmitry Voitsekhovskiy
 * All rights reserved.
 *
 * Redistribution and use in source and binary forms,
 * with or without modification, are permitted provided
 * that the following conditions are met:
 *
 * - Redistributions of source code must retain the above copyright notice,
 * this list of conditions and the following disclaimer.
 *
 * - Redistributions in binary form must reproduce the above copyright notice,
 * this list of conditions and the following disclaimer in the documentation
 * and/or other materials provided with the distribution.
 *
 * - Neither the name of the Moscow State University nor the names
 * of its contributors may be used to endorse or promote products
 * derived from this software without specific prior written permission.
 *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND
 * CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES,
 * INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
 * FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL
 * THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
 * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 ******************************************************************************/
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Linq.Expressions;

namespace Faml.Optimization
{
    [OptimizationContainer]
    public class UnfoldSubarrayOps
    {
        public static T[,,] UnfoldSubarray<T>(int width, int height, int depth, Expression<Func<int, int, int, T>> func, int[] ix,int[] iy,int[] iz)
        {
            T[,,] values = new T[ix.Length, iy.Length, iz.Length];
            Func<int, int, int, T> gen = func.Compile();
            for (int i = 0; i < ix.Length; i++)
                for (int j = 0; j < iy.Length; j++)
                    for(int k = 0;k< iz.Length;k++)
                        values[i, j, k] = gen(ix[i],iy[j],iz[k]);
            return values;
        }

        [OptimizationForPair("subarray", "unfold", 4, 0)]
        public static double[,,] UnfoldDoubleSubarray(int w, int h, int d, Expression<Func<int, int, int, double>> func, int[] ix, int[] iy, int[] iz)
        {
            return UnfoldSubarray<double>(w, h, d, func, ix, iy, iz);
        }
    }
}
